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ABSTRACT 

Wc  have  used  the  Infrared  Spectrograph  (IRS)  on  the  Spitzer  Space  Telescope  to  obtain  spectra  of  H  D  1 00764,  an 
apparently  single  carbon  star  with  a  circumstdlar  disk.  The  spectrum  shows  emission  features  from  polycyclic  aro¬ 
matic  hydrocarbons  (PAHs)  that  are  shifted  to  longer  wavelengths  than  normally  seen,  a  characteristic  of  “class  C” 
systems  in  the  classification  scheme  of  Peeters  et  a!.  All  seven  of  the  known  class  C  PAH  sources  are  illuminated  by 
radiation  fields  that  are  cooler  than  those  which  typically  excite  PAH  emission  features.  The  observed  wavelength 
shifts  are  consistent  with  hydrocarbon  mixtures  containing  both  aromatic  and  aliphatic  bonds.  We  propose  that  the 
class  C  PAH  spectra  arc  distinctive  because  the  carbonaceous  material  has  not  been  subjected  to  a  strong  ultraviolet 
radiation  field,  allowing  relatively  fragile  aliphatic  materials  to  survive. 

Subject  headings;  cireumstellar  matter  —  stars:  carbon 
Online  material:  color  figures 


I.  INTRODUCTION 

Complex  hydrocarbons  are  found  in  a  wide  variety  of  astro- 
physical  environments.  Gillcti  et  al.  (1973)  first  detected  what 
became  known  as  the  unidentified  infrared  (UIR)  emission  fea¬ 
tures,  the  strongest  of  which  are  at  3.3,  6.2,  7.7-7.9,  8.6,  1 1.3, 
and  12,7  pm.  Lcgcr  &  Puget  (1984)  and  Allamandola  et  al. 
( 1 985)  identified  the  carrier  of  these  features  as  polycyclic  aro¬ 
matic  hydrocarbons  { PAHs).  The  debate  over  PAHs  as  the  carrier 
of  the  UIR  features  raged  for  many  years,  and  while  the  identi¬ 
fication  now  seems  firm  (e.g.,  Allamandola  ei  al,  1999),  many 
questions  still  remain  over  the  origin,  evolution,  and  exact  com¬ 
position  of  astro  physical  hydrocarbons, 

Peeters  et  al.  (2002)  used  data  from  the  Short  Wavelength 
Spectrometer  (SWS)  on  the  Infrared  Space  Observatory  (ISO)  to 
group  PAH  spectra  into  three  classes  based  on  the  positions  of 
peak  emission  and  profiles  of  the  C-C  skeletal  deformation 
modes  within  the  PAHs  at  6.2  and  7. 7-7. 9  pm.  The  “class  A“ 
sources  display  features  at  6.22  and  7.65  pm  and  generally  in¬ 
clude  H  n  regions,  reflection  nebulae,  galaxies,  and  Herbig 
Ae/Be  ( HAeBe)  stars  still  embedded  in  the  regions  where  they 
formed.  The  “class  sources,  which  typically  include  plane¬ 
tary  nebulae  and  isolated  HAeBe  stars,  show  these  same  features, 
but  shifted  to  6,26  and  7.85  pm.  Peelers  et  al.  (2002)  also  iden¬ 
tified  two  “class  C”  sources,  including  the  well-known  Egg 
Nebula  (AFGL  2688),  which  are  in  transition  from  the  asymp¬ 
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totic  giant  branch  (AGB)  to  the  planetary  nebula  phase.  These 
two  sources  show  emission  at  6.26  jim,  no  emission  near  7.65  pm, 
and  a  broad  feature  centered  at  ^8. 2  pm.  It  remains  unclear  how 
these  different  classes  of  PAHs  arise  in  different  astrophysics] 
environments. 

Observations  with  the  Infrared  Spectrograph  ( IRS;  Houck  et  al. 
2004)  on  the  Spitzer  Space  Telescope  (Werner  et  al  2004)  have 
added  to  the  small  sample  of  class  C  PAH  spectra,  Kraemer  et  aL 
(2006)  reported  the  detection  of  a  class  C  PAH  spectrum  from 
MSX  SMC  029,  a  transition  object  in  the  Small  Magellanic  Cloud 
(SMC).  Jura  et  al.  (2006)  detected  a  class  C  PAH  spectrum  from 
HD  2335 1 7,  an  apparently  single  red  giant  (K2  III )  with  a  sub¬ 
stantial  infrared  excess,  most  likely  arising  from  a  circumstdlar 
disk  (Jura  2003).  This  system  is  remarkable  because  some  of  the 
hydrocarbons  were  probably  synthesized  from  the  oxygen-rich 
disk. 

HD  100764  is  one  of  a  handful  of  red  giants  (luminosity 
class  HI)  with  large  infrared  excesses,  characteristic  of  orbiting 
disks.  While  most  of  these  stars  are  oxygcn-rich  like  HD  233517, 
HD  100764  is  carbon-rich.  As  part  of  a  program  to  better  under¬ 
stand  these  stars,  we  obtained  IRS  data  for  HD  100764.  We  have 
found  that  its  IRS  spectrum  reveals  low-contrast  class  C  PAH 
features  not  detectable  in  previous  ground-based  spectroscopy. 
In  this  paper,  we  describe  the  PAH  spectrum  of  HD  100764  in 
detail  and  relate  it  to  other  class  C  PAH  sources. 

HD  100764,  also  known  as  CCS  1886  and  CGCS  3066,  is 
classified  as  an  early  carbon  star,  R2  (Sanford  1944)  or  Cl, I 
{ Yamashita  1972),  Using  optical  spectroscopy,  Dominy  ( 1984) 
found  that  the  central  star  has  a  temperature  of  4850  K,  is  iron 
poor,  and  has  an  enhanced  1 3C/  l2C  ratio,  making  this  carbon  star 
one  of  the  uncommon  J  stars.  The  Hipparcos  parallax  of2.78± 
1.18  mas  yields  an  uncertain  distance  of  ^360  pc  and  thus  a  total 
luminosity  of  73  Lv  (Parthasarathy  1991),  indicating  that  the  star 
is  a  first-ascent  red  giant. 

Lloyd  Evans  (1990)  identified  HD  100764  as  a  potential 
silicatc/earbon  star  because  of  its  red  [  12]  -  [25]  color,  as  measured 
by  the  Infrared  Astronomy  Satellite  (IRAS },  As  a  consequence,  it 
was  observed  twice  using  ground-based  mid-infrared  spectrome¬ 
ters  over  a  decade  ago.  LeVan  ct  al.  ( 1 992)  obtained  an  8- 1 4  pm 
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TABLE  1 

Low-Resolution  Wavelength  Ranges 


SL 

LL 

Order 

(/Jin) 

(/im) 

Second . . . 

.  5.10-7.59 

14.20-21.23 

Bonus . . . . 

.  7.73-8,39 

1 9.28— 2 1 ,23 

First  . 

.  7.59- 14.33 

20.46-37.00 

spectrum  at  the  Wyoming  Infrared  Observatory  (WIRO)  with  a 
resolution  (A/AA)  of  —50  and  a  signal /noise  (S/N)  ratio  of  —  5. 
While  they  could  rule  out  the  presence  of  silicate  dust  emission  at 
1 0  pm  in  the  spectrum,  they  could  say  little  more.  Skinner  ( 1 994), 
who  obtained  a  higher  quality  spectrum  with  CGS3  at  the  United 
Kingdom  Infrared  Telescope  (UKIRT),  also  saw  no  spectral  fea¬ 
tures  from  dust  and  concluded  that  the  infrared  emission  was 
dominated  by  amorphous  carbon  dust  in  a  circumstellar  disk.  The 
disk  explains  why  the  strong  infrared  excess  does  not  appreciably 
redden  the  optical  spectrum. 

Although  uncertain,  it  is  plausible  that  this  disk  formed  when 
HD  100764  became  a  red  giant  and  engulfed  a  companion.  The 
luminosity  of  HD  100764  is  consistent  with  other  R  stars  with 
measured  parallaxes,  which  Alksnis  ct  aL  (  1998)  have  deter¬ 
mined  are  first-ascent  red  giants,  not  AGB  stars.  McClure  (1997) 
observed  that  R  stars  do  not  have  companions  and  suggested 
that  their  carbon-rich  nature  and  the  lack  of  a  companion  could 
both  be  explained  by  coalescence  of  an  expanding  giant  and  its 
companion. 

2.  THE  INFRARED  SPECTRA 

2.1.  Observations  and  Analysis 

We  observed  HD  100764  with  both  the  Short-Low  (SL)  and 
Long-Low  (LL)  modules  of  the  IRS  in  standard  staring  mode  on 
2006  January  14.  Because  the  source  is  bright  for  the  IRS,  we 
used  the  shortest  possible  ramps  in  both  modules,  6  s  ramps  which 
provide  only  four  samples.  This  source  was  observed  indepen¬ 
dently  with  the  Short-High  (SH  )  and  Long-High  (LH)  modules 
in  Spitzer  program  278  (AOR  key  19481344),  and  we  include 
those  data  here  with  permission  of  the  PI,  M.  Wcmer. 

We  followed  the  standard  calibration  method  used  at  Cornell. 
We  started  with  the  S14  pipeline  output  from  the  Spitzer  Science 
Center  (SSC)  and  subtracted  a  sky  image  from  each  image  before 
extracting  a  spectrum.  The  high-resolution  observations  include 
separate  sky  images  obtained  dose  to  the  target.  For  SL,  the  sky 
image  for  one  spectral  aperture  is  obtained  when  the  source  is  in 
the  other  aperture  (i.c.,  aperture  differences  SL2-SL1  and  SL1  — 
SL2,  where  “SL2"  means  "SL  Order  2”).  For  LL,  we  simply 
differenced  the  images  with  the  source  in  the  same  aperture,  but 
in  the  other  nod  position.  For  a  source  as  bright  as  HD  100764, 
the  primary  benefit  of  background  removal  is  the  subtraction  of 
most  of  the  rogue  pixels.  These  pixels  have  temporarily  increased 
dark  currents  and  are  the  primary  source  of  noise  in  images  with 
bright  sources.  We  applied  imclean.pro7  to  replace  remaining 
rogue  and  bad  pixels. 

Spectra  were  extracted  from  the  individual  images  using  the 
SSC -defined  apertures  and  the  routines  profile,  ridge,  and 
extract,  which  are  available  with  the  SPICE  software  package 
(Spitzer/ IRS  Custom  Extraction).  To  calibrate  to  flux*density 
units,  we  used  a  spectral  correction  generated  from  IRS  obser- 

Available  from  the  SSC  as  irsclean  ■  pro. 
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Fig.  1. —  Spectrum  of  HD  100764  with  the  low-resolmion  modules  (black) 
and  with  the  high-rcsolution  modules  (gray).  [Sto*  the  electronic  edition  of  the 
Journal  for  a  color  version  of  this  figure.) 

vationsof  standard  stars:  HR  6348  ( K0  III )  for  SL,  and  HR  6348, 
HD  1 66780  (K4  III),  and  HD  173511  (K5  III  )  for  LL.  The  high- 
resolution  data  were  calibrated  using  £  Dra  (K2  III),  Error  bars 
are  based  on  the  standard  deviation  of  the  spectrum  in  the  two 
separate  nods  in  each  aperture.  The  resulting  spectrum  has  a  S/N 
ratio  typically  between  200  and  400. 

Wc  applied  scalar  multiplicative  corrections  to  the  segments  to 
remove  discontinuities  and  align  them  upwards  to  the  presum¬ 
ably  best-centered  spectral  segment,  in  this  case  LL2.  This  step 
required  corrections  of  9%  to  SL2,  2%  to  SLl ,  and  1%  to  LL1 , 
These  corrections  indicate  the  accuracy  of  the  pointing  of  the 
telescope  and  are  consistent  with  typical  observations.  Finally, 
wc  combined  the  bonus  order  with  valid  data  in  the  overlapping 
orders  and  removed  all  data  outside  the  ranges  of  trustworthy 
signal.  Table  1  gives  the  wavelength  ranges  used  for  each  low- 
resolution  order. 

The  flux  density  from  HD  100764  in  the  5-37  pm  spectral 
region  ranges  between  5  and  1 1  Jy,  and  we  have  therefore  as¬ 
sessed  the  saturation  and  nonlinearitics  in  the  response  function, 
especially  in  LL.  The  S14  pipeline  has  a  greatly  improved  cor¬ 
rection  for  these  problems,  and  it  produces  much  better  spectra 
compared  to  previous  pipeline  versions.  Figure  l  presents  the 
full  spectrum  of  HD  100764,  both  our  low-resolution  data  and 
also  the  high-resolution  data  obtained  by  M.  Werner  et  aL  A  com¬ 
parison  of  the  two  spectra  reveals  no  evidence  for  any  residual 
artifacts  from  the  saturation  and  linearity  corrections  in  our  low- 
resolution  spectrum.  One  artifact  that  does  remain  is  the  fringing 
from  the  LLI  filter  apparent  between  21  and  25  pm.  We  do  not 
attempt  a  correction,  because  the  analysis  in  this  paper  concen¬ 
trates  on  the  features  at  shorter  wavelengths.  For  this  paper,  we 
focus  on  the  low-resolution  spectrum  from  5  to  1 4  pm. 

Superimposed  on  the  continuum  produced  by  the  circumstel¬ 
lar  disk  are  emission  features  at  wavelengths  near  those  associ¬ 
ated  with  PAHs  and  a  broad  emission  feature  centered  at  10.6  ^m. 
To  isolate  the  PAH-related  features,  wc  subtracted  a  spline-fitted 
continuum,  as  Figure  2  shows.  This  spline  passes  through  and 
includes  the  broad  feature  peaking  at  1 0.63  /im.  We  examine  this 
feature  separately  in  §  2.5.  The  spectrum  of  HD  100764  also  shows 
a  narrow  absorption  band  from  acetylene  (C3H3)  at  13.7  pm. 
When  this  absorption  band  appears  with  PAH  emission  features, 
it  is  only  in  the  rare  class  C  PAH  spectra  { Kraemer  ct  al.  2006). 
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Fin*  2. —  Spectrum  of  HD  100764.  The  lop  panel  shows  the  co-added  spec¬ 
trum  and  the  spectra  from  the  individual  nod  positions  (offset  by  I  and  2  Jy,  re¬ 
spectively).  The  thin  solid  curve  is  a  spline  lit  to  the  data  at  wavelengths  marked 
with  diamonds.  These  points  lie  outside  of  the  PAH  features  andC^Hi  absorption 
band  at  1 3.7  /im.  The  bottom  panel  shows  the  resulting  continuum-subtracted  spec¬ 
trum,  along  with  the  error  bare  and  the  line  segments  fit  to  extract  the  strengths  and 
centers  of  the  spectroscopic  features.  The  line  segments  arc  thick  in  the  wavelength 
regions  used  to  lit  the  continuum  and  thin  where  the  feature  is  integrated.  The 
equivalent  width  of  the  13.7  acetylene  band  is  measured  before  any  contin¬ 
uum  subtraction.  {See  the  electronic  edition  of  the  Journal  fora  color  version  of 
this  figure.] 

The  low -contrast  features  in  the  8-14  pm  range  escaped 
previous  detection.  The  8  /mi  feature  has  a  contrast  of  5%,  but 
lies  at  the  edge  of  the  ground-based  window,  and  the  feature  near 
1 1 .3  /mi  has  a  contrast  of  only  3%,  too  small  to  be  seen  even  in 
the  spectrum  by  Skinner  ( 1 994),  which  has  a  S/N  ratio  of  ^  1 5  at 
those  wavelengths. 

Figure  2  (lop)  plots  the  spectra  from  the  two  nod  positions 
separately.  Systematic  errors  in  the  spectra  have  their  greatest 
impact  on  low-contrast  features,  and  they  often  show  up  as  dif¬ 


ferences  between  the  nods.  Some  detailed  differences  in  the 
PAH  features  are  apparent,  most  notably  at  ^8.5  pm,  but  these 
differences  arc  not  large  enough  to  have  a  significant  quanti¬ 
tative  effect  on  our  conclusions*  Figure  2  (bottom)  presents  the 
continuum-subtracted  spectrum  from  HD  100764,  along  with 
line  segments  used  to  extract  the  strength  and  position  of  each 
PAH  emission  feature,  as  described  in  §  2  J* 

2,2.  The  Comparison  Sample 

To  place  the  class  C  PAH  spectrum  from  HD  100764  in  con¬ 
text,  we  have  also  examined  several  other  PAH  spectra.  These 
comparison  spectra  include  the  two  class  C  PAH  sources  from 
the  SWS  database,  AFGL  2688  and  IRAS  13416-6243,  along 
with  SWS  spectra  of  NGC  1333  SVS  3  and  HD  44179  (the  Red 
Rectangle).  These  sources  are  prototypical  examples  of  class  A 
and  B  PAH  spectra,  respectively.  AH  four  spectra  are  from  the 
atlas  of  SWS  spectra  produced  by  Sloan  et  al.  (2003),  and  they 
have  been  regridded  to  a  lower  spectral  resolution,  in  the  case  of 
AFGL  2688,  the  spectrum  is  an  average  of  two  observations*  The 
comparison  sample  also  includes  the  four  sources  observed  by 
the  IRS  which  are  currently  known  to  have  PAH  spectra  either 
in  or  near  class  C,  including  the  K  giant  HD  233517  (Jura  et  al, 
2006),  the  post-AGB  source  MSX  SMC  029  (Kraemer  et  al. 
2006),  the  HAeBe  source  HD  135344  (Sloan  et  al.  2005),  and 
the  intenuedi ate- mass  T  Tauri  star  SU  Aur  (Furlan  et  al*  2006), 
Table  2  presents  the  comparison  sample,  including  their  spectral 
classes  and  effective  temperatures,  which  are  described  in  §  3.1 . 

2.3.  Extracting  the  PAH  Features 

For  each  source,  we  follow  the  algorithm  described  for  HD 
100764  by  spline  fitting  and  removing  a  continuum,  as  Figure  3 
illustrates.  It  has  proven  necessary  to  shift  the  spline  points  from 
one  spectrum  to  the  next  to  avoid  some  of  the  spectral  structure 
which  only  appears  in  individual  spectra.  Care  must  be  taken  in 
the  cases  of  IRAS  13416  and  AFGL  2688,  because  the  PAH 
features  have  such  low  contrast  with  respect  to  the  continuum. 
In  the  case  of  AFGL  2688,  the  features  at  1 2.0  and  1 2.7  /mi  must 
be  treated  with  caution,  because  their  apparent  strength  and  po¬ 
sition  shift  substantially  with  the  choice  of  spline  points.  The 
spline  fit  to  SU  Aur  removes  the  silicate  emission  feature;  this 
method  has  been  tested  for  several  similar  sources  showing  both 
PAH  and  silicate  emission  by  Keller  et  al*  (2007).  Two  of  the 
spectra  have  been  smoothed;  MSX  SMC  029  with  a  3  pixel  box¬ 
car  and  HD  2335 1 7  with  a  5  pixel  boxcar.  Figure  4  presents  the 
continuum-subtracted  PAH  spectra  for  HD  100764  and  the  com¬ 
parison  sample. 


TABLE  2 

HD  1 00764  anp  the  Comparison  Sample 


re¬ 
target  Instrument  A  OR  Key  or  TDT  Object  Type*  Spectral  Type*  (K)  PAH  Class 


HD  100764 . .  TRS  16262656  Red  giant  Cl,l  4850  C 

IRAS  13416-6243..*,.,..  SWS  62803904  Post-AGB  GI  I  5440  C 

AFGL  2688 . .....** .  SWS  33800604,35102563  Post-AGB  F5  la  6520  C 

MSX  SMC  029,,*™....  IRS  10656256  Post-AGB  ...  .  *  *  C 

HD  233517  . .. .  IRS  3586048  Red  giant  K2  III  4475  C 

SU  Aur....* . *,******.*..  IRS  3533824  T  Tauri  Gl  5945  C 

HD  135344.. . ..****..,.  IRS  3580672  Hcrbig  Ac/Be  F4  Ve  6590  B/C 

HD  44179*,* .  SWS  70201801  Post-AGB  A0  HI  9520  B 

NGC  1333  SVS  3....,-,,  SWS  65902719  Reflection  nebula  B6  14,000  A 


See  §  3. 1  for  details  and  references. 
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Fio.  3. —  Infrared  spectra  of  several  sources  showing  class  C  PAH  features, 
along  with  spectra  of  HD  44 1 79  (clais  B)  and  NGC  1333  SVS  3  (class  A),  The 
spectra  oflRAS  13416,  AFGL  2688,  NGC  1333  SVS  3,  and  HD  44 1 79  are  from 
the  ISOi SWS  database  (as  processed  by  Sloan  et  al  2003);  the  remaining  spectra 
arc  from  the  IRS.  The  thin  lines  show  the  con  tin  u  a  fit  to  the  spectra,  based  on  a 
spline  through  the  wavelengths  marked  with  diamonds.  The  resulting  continuum- 
subtracted  PAH  spectra  appear  in  Fig,  4  [See  (he  electronic  edition  of  the  Journal 
for  o  color  version  of  this  figure ,] 

Fitting  and  removing  a  spline  removes  most,  but  not  all,  of 
the  continuum  spectrum.  Because  we  have  not  forced  the  fitted 
spline  to  pass  precisely  through  the  continuum  on  either  side  of 
each  feature,  small  residual  tilts  remain.  While  these  would  have 
little  impact  on  the  extracted  strengths  of  the  features,  they  would 
shift  the  measured  center.  We  therefore  continue  to  use  the  method 
applied  to  a  sample  of  HAeBe  stars  showing  PAH  emission  by 
Sloan  et  al.  (2005).  As  they  did,  we  fit  line  segments  to  the  con¬ 
tinuum  on  either  side  of  each  PAH  feature  and  integrate  the  flux 
above  this  estimated  continuum  to  measure  the  strength  of  the 
feature.  We  determine  the  wavelength  with  half  the  emitted  flux 
to  cither  side,  which,  for  the  remainder  of  this  paper,  we  call  the 
central  wavelength. 

Other  authors  have  analyzed  the  position  of  the  peak  wave¬ 
length,  but  this  result  differs  from  the  central  wavelength  for  two 
reasons.  First,  PAH  features  are  usually  asymmetric,  so  the  peak 
will  lie  slightly  to  the  blue  of  the  center.  Second,  most  of  the 
features  are  blends  of  multiple  transitions.  The  feature  peaking  at 
^7. 7-7. 9  /xm  can  be  analyzed  as  a  combination  of  features  at 
7.65  and  7,85  /xm  (e.g,.  Peelers  et  al.  2002;  Bregman  &  Temi 
2005).  As  the  7,85  /xm  component  increases  with  respect  to  the 
7.65  /xm  component,  the  central  wavelength  will  move  steadily 
to  the  red,  but  the  peak  wavelength  will  not  move  significantly 
until  the  red  component  outshines  the  blue  component,  Analyz¬ 


Fia  4. — Continuum -subtracted  spectrum  of  HD  100764,  along  with  several 
comparison  spectra  (also  continuum -subtracted}.  The  name  of  each  source  is 
followed  by  the  classification  of  its  PAH  spectrum.  The  vertical  lines  are  at  the 
centers  of  the  PAH  features  for  a  class  B  spectrum;  6,26, 7.90, 1 1  26,  and  12,67  /xm. 
The  majority  of  the  class  C  spectra  show  features  shifted  to  the  red  of  these  lines. 
The  dashed  lines  at  6.85  and  7.25  mark  the  positions  of  aliphatic  emission  fea¬ 
tures.  The  bottom  two  spectra  are  the  prototypical  class  A  and  class  B  spectra 
from  the  /SO/SWS  database.  The  remainder  are  the  sample  of  known  class  C  PAH 
spectra  (or  in  the  case  of  HD  135344,  nearly  class  C  PAH  spectrum  ).  The  extrac¬ 
tions  oflhe  PAH  features  at  1 2,0  and  1 2.7  jim  in  AFGL  2688  arc  particularly  un¬ 
trustworthy  due  to  Lhcir  extremely  low  contrast  and  resultant  sensitivity  to  the 
location  of  the  spline  points  for  the  continuum.  In  the  spectrum  of  AFGL  2688,  the 
bottom  of  the  1 3 ,7  feature  has  been  truncated  to  prevent  it  from  overrunning 
the  spectra  below,  [See  the  electronic  edition  of  the  Journal for  a  color  version  of 
th  is  figure,] 

ing  the  central  wavelength  provides  better  insight  about  the  rel¬ 
ative  strengths  of  the  various  components  of  each  PAH  feature. 

Table  3  gives  the  wavelength  ranges  used  to  fit  the  line  seg¬ 
ments  to  each  feature;  these  have  changed  some  from  the  ranges 
used  by  Sloan  et  ah  (2005).  We  use  the  class  C  wavelengths  for 
all  of  the  sources  except  HD  441 79  and  NGC  1 333  SVS  3,  which 
require  shifted  wavelengths  to  extract  the  features.  All  of  the  spec¬ 
tra  except  HD  100764  and  SU  Aur  have  a  feature  at  8.6  /xm  from 
the  C-H  in- plane  bending  mode  on  the  red  shoulder  of  the  C-C 
feature  peaking  between  7.7  and  8.2  /xm.  To  eliminate  the  influ¬ 
ence  of  the  8,6  /xm  feature  on  the  position  and  strength  of  the 
broader  feature,  we  fit  a  line  segment  underneath  it  and  include 
only  the  flux  below  that  line  in  our  measurement  of  the  7,7-  8,2  /.xm 
feature. 

Table  4  presents  the  extracted  central  wavelengths,  while 
Table  5  gives  the  equivalent  fluxes  (i,c.3  integrated  flux  densities), 
and  for  the  acetylene  band  at  13.7  /xm,  the  equivalent  width.  In 
cases  where  no  feature  is  apparent  or  the  strength  of  the  extracted 
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TABLE  3 

Wavelength  Intervals  Used  to  Extract  Spectral  Features 

For  Class  C  Sources  For  Class  A  and  B  Sources 


Feature  ^biue  ^«d  ^™3 

(F  m)  (ftm)  (^m)  (jjm)  U*m) 


6.2 . * . ,*  5.80-5.95  6.65-6.80  As  class  C  6*54-6.69 

7.9.. ....... . *******  7.30-7.60  9.00-930  7  10-7,40  8.92-9.17 

B.6..,**** . * . **  832-8  44  8.86-8,98  As  class  C  As  class  C 

113 . . . . .  10.80-11.00  H. 85- 11.95  10.68-10.86  11*72-11*90 

12.0 . „* . . . ********  11,85-11*95  1235-12*60  11*72-11.90  12.13-12*25 

12.7.. ** . . . . .  12.35-12.6Q  13.25-13,40  12,07-12*19  12,98-13,10 

13.7 .  13.25-13.55  13.83-14.13  As  class  C  As  class  C 


feature  has  a  S/N  ratio  less  than  l ,  Table  4  does  not  quote  a  central 
wavelength.  Table  4  includes  notes  indicating  those  wavelengths 
for  features  with  S/N  ratios  less  than  2  or  3.  Results  for  some  of  the 
comparison  spectra  have  been  published  recently,  but  the  results 
presented  here  represent  newer  S 14  pipeline  output  and  reflect  mod¬ 
ifications  to  the  continuum  fitting  and  extraction  wavelengths. 

2,4*  The  PAH  Spectra 

Figures  3  and  4  include  the  spectra  of  our  class  A  and  B  pro¬ 
totypes,  NGC  1333  SVS  3  and  HD  441 79.  Using  these  and  other 
SWS  data,  we  find  that  the  C-C  modes  produce  features  in  class 
A  spectra  with  typical  central  wavelengths  of  6.23  and  7*7  pm, 
while  class  B  spectra  have  features  centered  at  6,26  and  7.9  /*m 
(in  essential  agreement  with  Pecters  ct  al.  2002).  In  the  two  class 
C  spectra  in  the  SWS  database,  these  features  are  shifted  to  6.30 
and  8.2  pm,  respectively.  In  HD  100764,  they  arc  centered  at 
6.33  and  8, 1 5  /an.  In  the  other  class  C  IRS  sources,  these  features 
appear  in  the  6,25-6.28  ft m  and  8.09-8.29  ftm  ranges. 

The  features  beyond  10  pm  arise  from  out-of-plane  C-H 
bending  modes.  The  wavelength  of  these  modes  depends  on  the 
number  of  adjacent  hydrogen  atoms  on  a  single  aromatic  ring 
(e.g*,  Atlamandola  et  al.  1 989).  Single  H  atoms  produce  tile  solo 
mode,  which  class  A  and  B  spectra  show  between  I  1.21  and 
1 1*30  /tm,  with  a  mean  at  1  1.26  fi m*  Rings  with  two  adjacent 
H  atoms  produce  the  duo  mode,  which  normally  appears  at 
12.00  ;j.m.  Three  adjacent  H  atoms  produce  the  trio  mode,  which 
appears  at  1 2,70  pm  in  the  class  A  and  B  spectra. 

In  alt  of  the  class  C  sources,  the  solo  out-of-plane  bending 
mode  has  shifted  to  longer  wavelengths.  In  the  comparison  class 
C  PAH  spectra  observed  by  the  IRS,  the  features  have  shifted  from 
the  nominal  1 1 *26  ftm  position  to  between  1 1 .32  and  1 1.36  /im. 


In  the  two  class  C  spectra  in  the  SWS  database,  the  shifts  are  dif¬ 
ficult  to  measure,  due  to  possible  effects  from  molecular  absorp¬ 
tion  and  artifacts*  Our  best  estimate  places  the  band  centers  in  the 
11.36-1 1.42  pm  range.  In  HD  100764,  this  feature  is  centered 
further  to  the  red,  at  1 1 ,47  fim. 

The  duo  and  trio  modes  also  appear  to  shift  in  the  class  C  spec¬ 
tra,  but  they  are  fainter  and  more  difficult  to  measure.  A  shift  is 
apparent  in  both  features  in  the  SWS  spectra  of  AFGL  2688  and 
IRAS  13416,  from  12,0  to  — 12.1  — 12.2  /im  for  the  duo  mode  and 
from  12.7  to  M2.9  /on  for  the  trio  mode,  but  these  results  should 
be  treated  with  reservation.  The  problem  is  that  the  contrast  of 
these  features  is  small,  and  their  apparent  strength  and  position 
depend  on  the  details  of  how  we  fit  a  spline  to  the  continuum. 
Of  the  IRS  comparison  spectra,  SU  Aur  shows  the  dearest  duo¬ 
mode  feature,  but  the  S/N  ratio  is  only  —2  and  the  shift  is  minimal. 
The  spectrum  of  HD  1 00764  shows  a  more  significant  shift,  but 
again,  the  S/N  ratio  is  only  —2*  The  situation  with  the  trio  mode 
is  similar  in  the  IRS  spectra,  A  shift  is  apparent,  but  the  uncer¬ 
tainties  in  central  wavelength  arc  large. 

In  summary,  the  positions  of  the  PAH  features  in  the  spectrum 
of  HD  100764  are  fully  consistent  with  a  class  C  PAH  spectrum. 
Two  of  the  features,  at  6*33  and  1 1 ,47  pm,  represent  extreme  po¬ 
sitions  for  our  sample.  The  spectrum  is  definitely  not  class  A  or  B. 

2*5,  Other  Spectral  Features 

In  Figures  2  and  4,  the  spectrum  of  HD  100764  shows  struc¬ 
ture  between  the  PAH  features  at  6,3  and  8.1  fim  reminiscent  of 
the  features  observed  in  emission  at  6.85  and  7.25  ftm  in  the 
spectra  of  the  HAeBe  stars  HD  34282  and  HD  1 69142  by  Sloan 
et  al.  (2005).  These  features  are  most  obvious  in  the  spectrum  of 
AFGL  2688  (where  they  can  be  seen  in  both  scan  directions). 


TABLE  4 

Central  Wavelengths  of  the  PAH  and  QHj  Features 


L  (/i m) 


Tahoet  6.2  /an  PAH  7.7-K.2  /<m  PAH  8.6  PAH  ]  1.3  /im  PAH  12.0  ,im  PAH  12.7  ,iro  PAH  13.7  ,tm 


HD  100764 .  6.33  ±  0.01  8,15  ±  0.06  ...  11.47  ±  0.07  12.14  ±0.15*  12.83  ±  0.27*  13,70  ±0.02 

[R AS  13416 . 6.31  ±0.01  8,22  ±  0.02  8.71  ±  0.02  11.42  ±  0.03  12.25  ±  0,07  12.87  ±0.10  13.70  ±0.00 

AFGL  2688 .  6.29  ±  0.01  8.20  ±  0.0 1  8.68  ±  0.03  11.36  ±  0.03  12.13  ±  0.03  12.93  ±  0.03  13.71  ±  0.01 

MSXSMC  029 .  6.25  ±  0.04  8.12  ±  0.10  8,68  ±  0.05  11.32  ±  0.07  ...  13.67  ±  0,02 

HD  233517 .  6.28  ±  0.03  8.26  ±  0.03  8,72  ±  0.07  11.35  ±  0.02  ...  12.75  ±  0.09 

SU  Aur .  6.26  ±0.01  8.19  ±  0.03  8.67  ±0.18*  11.33  ±  0.0 1  12.02  ±  0.03b  12.78  ±  0.08 

HD  135344 .  6,28  ±  0.03  8.08  ±  0.04  8.64  ±  0.07  11,35  ±  0.01  ...  12.75  ±  0.08 

HD  44179 .  6.26  ±  0.00  7.84  ±  0,06  8,66  ±  0,01  11.27  ±  0.00  12.00  ±  0.04  12.71  ±  0.02 

NGC  1333  SVS  3 .  6.23  ±  0.00  7,68  ±  0.01  8,61  ±  0.01  11.26  ±  0.01  12.00  ±  0.02  12.70  ±0.01 


#  S/N  of  extracted  feature  <2, 
b  S/N  of  extracted  feature  <3, 
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TABLE  5 

Extracted  Strengths  of  the  PAH  and  C2H2  Features 

F(10-'5  Wm2)  Wi(\0r*itmr 


Targct  6.2  pm  7.T-8.2  pm  8.6  ftm  1 1.3  pm  12.0  pm  12.7  pm  B.7  pm 


HD  100764 . . . . .  4.61  ±  0.24  9.15  ±  0,82  „  .  2.27  ±  0.28  0.32  ±0.17  0.16  ±  0,13  3.18  ±  0.64 

IRAS  13416. . 57.3  ±  2 A  225  ±5  1.73  ±  0.51  22.9  ±  1.5  7.39  ±  L66  7.71  ±  1.54  5.10  ±  0.48 

AFGL  2688.. . MI  ±  4  860  ±  10  18.6  ±  1.6  147  ±  15  76.5  ±  7,4  204  ±  14  6,23  ±  0.57 

MSX  SMC  029 . . .  0. 16  ±  0,03  0,29  ±  0.03  0.051  ±  0.007  0.068  ±  0.012  0.005  ±  0,021  0.001  ±0,015  14.84  ±  2.96 

HD  233517 . . . .  1,46  ±0,15  1.75  ±0.08  O.Jt  ±  0.03  0.66  ±  0.04  ...  0.091  ±  0,037  1.22  ±  3.09 

SU  Aur . ,i . . .. .  6,08  ±  0.35  6.64  ±  0.35  0.14  ±  0.12  3.77  ±  0.13  0,080  ±  0,037  0.31  ±  0.08 

HD  135344 . , . . . .  5.44  ±  0.54  9,42  ±  0.46  0.43  ±  0J2  2.34  ±  0.06  ...  0. 14  ±  0.04 

HD  44179  . . 2969  ±  23  5731  ±  677  960  ±19  1182  ±  19  48,3  ±  10.6  239  ±  14 

NGC  1333  SVS  3 . 241  ±  3  280  ±  7  69.9  ±1.8  172  ±  5  3.52  ±  0.63  57,2  ±  1,2 


"  Equivalent  width  of  the  acetylene  absorption  band. 


The  carriers  of  the  6,85  and  7.25  pm  features  have  been  iden- 
lifted  as  aliphatic  hydrocarbons  (Chiar  et  ah  2000)  and/or  hy¬ 
drogenated  amorphous  carbon  ( HAC;  Furton  et  ah  1 999).  In  HD 
1 00764,  the  features  are  faint  and  close  to  our  uncertainty  level. 
The  7.25  /Jtn  feature  is  stronger,  which  is  a  little  unusual.  As  a 
consequence,  our  identification  of  these  features  in  HD  100764 
is  only  tentative.  Other  class  C  spectra  also  show  possible  spec¬ 
tral  structure  in  this  area,  but  with  some  variations.  For  example, 
IRAS  13416  shows  an  extra  feature  {all  three  appear  in  both  scan 
directions).  The  class  A  and  B  prototypes  do  not  show  the  6.85 
and  7.25  pm  features. 

Excluding  the  emission  features  from  PAHs  and  related  mol¬ 
ecules,  the  remaining  spectrum  from  HD  100764  consists  of  a 
cool  dust  continuum,  absorption  due  to  acetylene  at  13.7  /tm  and 
a  broad  emission  feature  centered  at  M0.6  pm.  The  acetylene 
band  is  centered  at  13.7  /on,  which  is  typical  for  a  carbon  star  in 
the  Galactic  sample,  and  it  has  an  equivalent  width  of  only 
0.0032  ±  0.0006  ^m,  about  a  factor  of  10-20  times  weaker  than 
seen  in  the  SWS  sample  of  Galactic  carbon  stars  (Sloan  ct  al. 
2006).  The  acetylene  bands  arc  apparent  in  sources  with  opti¬ 
cally  thick  dust  shells,  indicating  that  the  absorbing  molecules 
arc  not  in  the  stellar  photosphere,  but  instead  intermixed  w  ith  or 
above  the  emitting  dust  (Sloan  et  al.  2006  and  references  within). 
In  the  case  of  HD  100764,  it  is  similarly  likely  that  the  acetylene 
absorption  is  not  from  the  stellar  photosphere.  Instead,  the  acety¬ 
lene  is  probably  distributed  in  and/or  above  the  disk. 

The  broad  emission  feature  to  the  blue  of  the  1 1 .47  pn\  PAH 
feature  can  be  fit  reasonably  well  with  a  Gaussian  centered  al 
10.63  jim  with  a  full  width  at  half  maximum  of  1 .37  pm.  As  Speck 
et  al.  (2005)  have  explained,  typical  SiC  dust  grains  emit  at 
^  1 1.3  pm*  but  larger  SiC  grains  emit  predominantly  at  1 0.8  pm^ 
because  the  contribution  from  the  longitudinal  optic  at  10,8  pm 
increases  with  respect  to  the  transverse  optic  at  II. 5  pm  as  the 
grain  size  grows,  Sloan  et  al  ( 2006}  may  have  observed  this  effect 
in  two  carbon  stars  in  the  SMC,  but  they  were  trying  to  explain  a 
shift  from  1 1 .3  to  1 1 .0  pm.  A  shift  past  1 0,8  //m  to  1 0,6  pm  seems 
far  less  likely.  We  conclude  that  this  broad  emission  feature  at 
10.6  pm  is  not  likely  to  arise  from  SiC;  it  remains  unidentified. 

3.  DISCUSSION 

3.L  Dependence  of  the  PAH  Emission 
on  Excitation  Temperature 

Table  2  gives  the  spectral  types  and  effective  temperatures  of 
our  PAH  sample.  HD  100764  has  been  classified  as  R2  or  C1J 
(Sanford  1944;  Yamashtta  1972)  and  has  an  effective  tempera¬ 
ture  of  4850  K  { Dominy  1 984),  The  two  class  C  sources  in  the 


SWS  database  appear  to  be  post-AGB  objects  evolving  into  plan¬ 
etary  nebulae  (PNe),  The  central  star  in  AFGL  2688  is  classified 
as  F5  la,  and  the  surrounding  nebulosity  shows  no  strong  emis¬ 
sion  lines  (Crampton  et  al,  1975).  IRAS  1341 6-6243  has  a  spec¬ 
tral  type  of  G1  1  (Hu  et  al.  1993)  and  also  shows  no  ionized 
emission  ( Kaufl  etaL  1 993). 8  Van  der  Veen  et  al.  (1989)  describe 
it  as  in  transition  from  the  AGB  to  a  PN.  We  have  determined 
temperatures  from  the  spectral  classes  for  these  two  objects  using 
the  relation  calibrated  by  Me  William  (1991). 

The  PAH  spectrum  of  MSX  SMC  029  is  best  described  as 
class  C.  While  the  6.2  /im  feature  has  not  shifted  much  from  the 
position  seen  in  class  B,  the  7-9  pm  PAH  complex  and  11.3  pm 
feature  are  shifted  to  die  red.  Like  the  post-AGB  objects  AFGL 
2688  and  IRAS  13416,  MSX  SMC  029  shows  an  absorption 
band  from  acetylene  at  1 3.7  pm.  Kraemer  et  al.  (2006)  argue  that 
it  is  one  of  the  few  known  post-AGB  objects  in  the  SMC.  It  is 
unfortunate  that  there  are  no  optical  data  to  help  classify  the 
spectrum  of  the  central  source. 

Sloan  et  al.  (2005)  presented  the  spectra  of  four  HAcBc  stars, 
three  classified  as  early  to  mid  A  dwarfs,  and  one,  HD  135344, 
as  F4  Ve  (Dunkin  ct  al.  1997),  which  corresponds  to  an  effec¬ 
tive  temperature  of  6590  K  (using  the  calibration  of  Kenyon  & 
Hartmann  1995,  Table  A5),  The  PAHs  around  the  A  dwarfs  re¬ 
semble  class  B  spectra,  except  that  the  7.85  pm  feature  is  shifted 
to  7. 9-8,0  /im,  and  in  the  F  dwarf  HD  135344,  this  feature  is 
shifted  even  further,  to  8.08  pmy  making  it  intermediate  between 
class  B  and  C.  Keller  et  al.  (2007)  have  re-examined  these  four 
spectra  as  par!  of  a  larger  sample  of  1 8  HAeBe  stars  and  four 
related  stars,  and  they  confirm  these  previous  results.  As  a  group, 
the  HAeBe  stars  exhibit  PAH  spectra  which  differ  from  the  typ¬ 
ical  class  B  PAH  spectra,  with  more  redshifted  features.  The 
coolest  stars  in  the  sample  have  spectra  which  are  even  closer 
to  class  C.  One  source  examined  by  Keller  ct  al.  (2007)  is  the 
intermediate-mass  T  Tauri  star  SU  Aur,  which  has  a  spectral  class 
of  G 1  and  an  effective  temperature  of  5945  K  (Calvet  ct  ah  2004), 
Furlan  et  at.  (2006)  noted  that  while  silicate  emission  dominates 
its  infrared  spectrum,  emission  from  PAHs  is  apparent  at  6.2  pm 
and  on  the  wings  of  the  10  pm  silicate  feature  at  8  and  1 1 .3  pm. 
This  source  is  the  coolest  star  in  the  sample  of  Keller  ct  al.  (2007), 
and  it  is  their  one  source  with  an  unambiguous  class  C  PAH 
spectrum  L 

HD  2335 17  has  a  spectral  type  of  K2  III  (  Fekdctal.  1996),  an 
effective  temperature  of 4475  K  (Balachandran  et  al.  2000),  and 
in  the  infrared,  a  class  C  PAH  spectrum.  It  is  the  coolest  object  in 

*  Van  de  Sieenc  cl  al.  ( 2000)  actually  detected  Bn  absorption  al  2. 17  pm. 
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Fig.  5, — Correlation  between  the  centra]  wavelength  (defined  m  §  2,3)  of  the 
two  strongest  PAH  features  and  the  effective  temperature  of  the  exciting  star  The 
top  panel  plots  the  position  of  the  7J-8.2  pm  PAH  complex,  and  the  bottom 
panel  plots  the  position  of  the  11 .3  ftm  feature.  This  figure  contains  three  Hcrbig 
Ae/Bc  stars  not  included  in  Tables  2, 4,  or  5,  The  IRS  data  for  these  sources,  HD 
34282,  HD  14 1 569,  and  HD  169142,  appear  elsewhere  (Sloan  el  al.  2005;  Keller 
cl  al.  2007).  HD  1 35344  is  the  right-most  square.  The  two  prototypes  are  NGC 
1333  SVS  ^{  far  left)  and  HD  44 1 79  {second  to  left).  [See  the  electronic  edition  of 
the  Journal  for  a  color  version  of  this  figure.] 


our  comparison  sample,  and  the  center  of  its  7-9  pm  PAH  com¬ 
plex  is  also  at  the  longest  wavelength  observed. 

By  comparison,  the  class  A  and  B  prototypes  are  warmer,  NGC 
1 333  SVS  3  has  a  spectral  class  of  B6  (Harvey  et  aE.  1 984),  which 
corresponds  to  an  effective  temperature  of  14,000  K  (Kenyon  & 
Hartmann  1995).  HD  44179  has  a  spectral  type  of  A0  III  (Cohen 
eta!.  197  5),  corresponding  to  9520  K(Kcny  on  &  Hartmann  1995). 
While  a  range  of  temperatures  have  been  used  to  model  the  Red 
Rectangle,  we  adopt  9520  K  as  most  consistent  with  the  optical 
spectral  class. 

Figure  5  plots  the  positions  of  the  PAH  features  as  a  function 
of  the  effective  temperature  of  the  exciting  source.  In  addition  to 
the  comparison  sample  in  Tables  2,  4,  and  5,  the  figure  also  in¬ 
cludes  the  HAeBe  stars  HD  34282,  HD  141569,  and  HD  169142. 
Sloan  ct  al.  ( 2005)  first  published  the  ERS  spectra  of  these  sources, 
and  the  results  presented  here  are  from  the  more  recent  analysis 
by  Keller  et  af  (2007),  The  correlation  of  the  position  of  the 
7.7 -8, 2  /mi  PAH  complex  with  effective  temperature  is  impres¬ 
sive,  While  the  1 1 .3  /mi  PAH  feature  also  tends  to  shift  to  longer 
wavelength  as  the  temperature  decreases,  the  trend  is  not  so  clear, 
with  little  difference  between  the  class  A  and  B  sources  and  more 
scatter  in  the  class  C  sources. 

3.2.  Spatial  Dependence  of  the  PAHs  in  Extended  Sources 

In  the  original  PAH  model,  ultraviolet  radiation  was  required 
to  excite  the  U1R  emission  features  (e.g.,  Allamandola  et  al. 
1989).  Uchidaet  al,  ( 1998,  2000)  first  observed  PAH  emission 
in  cooler  radiation  fields,  with  effective  temperatures  as  low  as 


6800  K  in  the  case  of  the  reflection  nebula  vdB  133,  While  these 
observations  initially  cast  some  doubt  on  the  PAH  model,  Li  <& 
Draine  (2002)  showed  that  the  absorption  of  visible  photons  by 
large  or  ionized  PAHs  can  produce  the  observed  UIR  features. 

Since  PAH  emission  often  arises  in  a  photodissociation  region 
(PDR),  where  a  front  is  eating  its  way  into  a  denser  region  and 
destroy ing  the  PAHs  in  the  process,  one  might  expect  a  change  in 
the  nature  of  the  PAH  spectrum  as  the  temperature  of  the  exciting 
radiation  field  varies.  Uchida  el  al,  (2000)  did  not  observe  varia¬ 
tion  in  the  relative  strengths  of  the  features  in  their  sample,  but  they 
did  observe  a  change  in  the  shape  of  the  7-9  pm  PAH  complex. 

Bregman  &  Tend  (2005)  analyzed  the  shape  of  the  7-9  pm 
complex  as  a  function  of  position  within  the  reflection  nebulae 
vdB  133,  NGC  1333  SVS  3,  and  vdB  1 7  (also  within  NGC  1333). 
They  found  that  closer  to  the  exciting  star,  the  strength  of  the 
7.65  pm  component  increases  relative  to  the  7.85  pm  compo¬ 
nent,  thus  shifting  the  measured  center  of  the  7.7-7.9  pm  PAH 
feature  to  the  blue.  They  also  found  that  this  complex  grows 
stronger  relative  to  the  C-H  modes  at  longer  wavelengths  near 
the  exciting  sources  and  suggested  that  the  wavelength  shift  at 
7.7-7.9  pm  might  be  related  either  to  the  ionization  fraction  of 
the  PAHs  or  to  PAH  processing  by  the  radiation  field,  as  previ¬ 
ously  proposed  by  Hony  ct  al,  (200 1 ), 

3.3.  The  Possible  Rote  of  Aliphatics 

Hydrocarbons  can  exist  in  both  aromatie  or  aliphatic  forms. 
Aromatic  hydrocarbons  are  dominated  by  sp2  bonds  which  pro¬ 
duce  the  polycyclic  ring  structure  that  define  PAHs.  Aliphatics 
include  sp 1  bonds,  as  seen  in  acetylene,  and  sp*  bonds,  as  seen  in 
methane  and  longer  alkane  chains.  The  sp2  bonds  efficiently 
spread  the  energy  absorbed  from  a  photon  over  the  entire  mol¬ 
ecule  before  it  can  disrupt  a  single  bond,  making  the  aromatic 
carbon  structure  in  PAHs  more  stable  than  aliphatic  bonding 
structures.  Thus,  when  a  mixture  of  aliphatic  and  amorphous  hy¬ 
drocarbons  is  exposed  to  a  harsh  radiation  environment,  we  ex¬ 
pect  the  aliphatic  bonds  to  be  broken  first.  We  propose  that  the 
class  C  spectra  arise  from  mixtures  of  aromatic  and  aliphatic  ma¬ 
terial  that  have  not  yet  been  subjected  to  intense  ultraviolet  radia¬ 
tion,  allowing  the  more  fragile  aliphatic  bonds  to  survive. 

Consider  the  class  C  spectrum  of  MSX  SMC  029,  This  source 
also  shows  a  strong  absorption  band  from  acetylene  (C2H2)  at 

13.7  /an,  centered  within  a  broader  absorption  band  in  the  12- 
16  /*m  region  (see  Kraemcr  et  al.  2006,  Fig.  4).  Kracmer  et  al. 
(2006)  argue  that  this  broader  band  is  due  to  other  aliphatic  mol¬ 
ecules  seen  in  a  similar  spectrum  at  higher  spectral  resolution  in 
SMP  LMC  1 1,  a  source  in  the  Large  Magellanic  Cloud  in  tran¬ 
sition  to  the  early  planetary  nebula  stage  (Bemard-Salas  et  al 
2006,  Fig.  1).  Acetylene  appears  in  other  class  C  PAH  spectra, 
including  AFGL  2688,  IRAS  13146,  and  our  spectrum  of  HD 
10 07 64. 9  The  association  of  the  acetylene  absorption  band  at 

1 3 .7  pm  and  the  broader  12-16  pm  absorption  with  elass  C  PAH 
emission  is  consistent  with  our  proposal  that  the  class  C  spectra 
arise  from  mixtures  of  PAHs  and  aliphatics. 

The  ground-based  mid-infrared  spectrum  of  HD  56126  (IRAS 
07134+1005)  shows  rcdshifted  PAH  features  similar  to  our  IRS 
spectrum  of  HD  100764,  with  emission  both  in  a  plateau  at 
8.0-8.2  /mn  and  from  out-of-plane  bending  modes  at  1 1.4  and 
12.2  pm  (iusttanont  et  al,  1996),  Although  not  as  pronounced, 
the  infrared  spectra  of  IRAS  04296+3429  and  IRAS  05341  + 
0852  display  similar  unusual  characteristics.  All  three  ofthese 
sources  are  post-AGB  objects.  HD  56126  is  classified  as  F5  I 


y  It  is  missing  from  HD  2335 1 7,  which  is  an  oxygen-rich  source,  and  the  disks 
around  HD  135344  and  SU  Aur. 
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(Nassau  et  al.  1965),  while  IRAS  05341  and  IRAS  04296  areF4 
lab:  and  GO  la,  respectively.  These  sources  also  show  unusual 
spectra  in  the  3  ft m  range,  which  is  dominated  by  C-H  stretching 
modes.  IRAS  05341  shows  a  feature  at  3*40  jim  that  dominates 
the  normally  stronger  feature  at  3.29  j*m  (Geballe  &  van  der 
Veen  1990),  while  IRAS  04296  also  shows  an  unusually  strong 
3.40  fim  component  (Geballe  ct  al,  1992).  The  3*29  ftm  feature 
arises  from  an  aromatic  C-H  stretching  mode  (e.g.,  Duley  & 
Williams  1981;  Allaniandola  ct  al.  1989),  while  the  3.4  fi m  fea¬ 
ture  arises  from  a  similar  mode  in  aliphatic  C-H  bonds  ( Duley  & 
Williams  1981;  Jourdain  de  Muizon  et  al  1990;  Geballe  et  ah 
1994). 

Geballe  et  ah  (1989)  observed  an  increase  in  the  ratio  of  the 
aliphatic  3.4  fim  feature  to  the  aromatic  3.29  /xm  feature  in  the 
shielded  molecular  zone  of  the  Orion  Bar  compared  with  the  ion¬ 
ized  zone.  Sloan  et  al.  ( 1997)  verified  these  observations  with 
long-slit  3  /xm  spectroscopy  and  showed  that  the  3.4  /xm  feature 
actually  arises  from  multiple  aliphatic  components.  Joblin  et  aL 
( 1996)  obtained  spatially  resolved  3  /mi  spectroscopy  in  the  re¬ 
flection  nebulae  NGC  1 333  and  NGC  2023  and  found  that  where 
the  ultraviolet  radiation  field  was  less  intense,  the  3.4  ftm  feature 
was  stronger  compared  to  the  3.29  pm  feature.  All  of  these  authors 
argued  that  the  radiation  field  is  altering  an  initial  mix  of  ali¬ 
phatic  and  aromatic  bonds  by  destroying  the  more  fragile  ali¬ 
phatic  bonds. 

3.4.  Hydrogenated  Amorphous  Carbon 

Hydrogenated  amorphous  carbon  (HAC)  is  a  generic  name 
for  a  mixture  of  aliphatic  and  aromatic  carbon,  consisting  of 
PAH  dusters  embedded  within  a  matrix  of  aliphatieally  bonded 
material  (see  Jones  et  al.  1990  and  references  therein)*  This  de¬ 
scription  is  similar  to  the  material  we  arc  proposing  to  explain  the 
class  C  PAH  spectra. 

Compared  to  previously  obtained  class  C  PAH  spectra,  HD 
100764  represents  an  extreme  case,  because  its  features  are  more 
consistently  rcdshiflcd  and  its  exciting  radiation  field  is  one  of 
the  coolest  known.  The  PAHs  around  this  source  must  be  rela¬ 
tively  unprocessed  by  energetic  photons,  making  this  source  an 
excellent  candidate  for  the  presence  of  HAC-like  structures. 

Figure  6  shows  two  representative  laboratory  spectra  of  a 
HAC  sample  at  ^800  K  (Scott  1997).  HAC  films  with  a  thick¬ 
ness  of~l  jim  were  prepared  by  ablating  graphite  in  a  hydrogen 
atmosphere  with  an  cxcimer  laser,  as  reported  previously  (Scott 
&  Duley  1996;  Scott  et  al.  1997a)*  Films  were  deposited  on  a 
stainless  steel  substrate,  placed  in  a  vacuum  chamber,  and  heated 
by  conduction  to  ^800  K.  Infrared  emission  spectra  were  ob¬ 
tained  with  a  Fourier-transform  spectrometer  at  a  resolution  of 
4  cm  1  over  a  spectral  range  of  500-6000  cm"1  (1.7-20  pm). 
The  spectra  were  calibrated  by  taking  the  ratio  of  the  emission 
from  H  AC  to  that  of  an  uncoated  stainless  steel  substrate  at  the 
same  temperatures  and  making  a  cubic  polynomial  baseline  cor* 
rection.  The  two  spectral  segments  from  5.5  to  8.0  pm  and  7.5  to 
14  pm  in  Figure  6  were  obtained  separately,  and  their  relative 
intensities  arc  unknown.  We  have  simply  scaled  both  segments 
to  a  maximum  of  l  .0  arbitrarily  in  the  figure. 

While  the  relative  band  strengths  in  Figure  6  are  not  well 
determined,  the  positions  are  more  clear.  The  emission  features 
have  shifted  to  the  red,  with  the  nominal  6.23  and  IL26  pm 
features  appearing  at  6.32  and  1 1 .45  jirn,  respectively.  If  the  fea¬ 
tures  in  the  12-13  pm  range  can  be  attributed  to  the  duo  and 
trio  C-H  out-of-plane  bending  modes,  then  they  have  shifted 
from  1 1.99  and  12.67  /mi  to  12*39  and  13.43  /tm.  The  feature 
centered  at  8.53  /mi  could  be  a  red  shifted  C-C  mode,  or  it  could 
simply  be  the  C-H  in-planc  bending  mode.  This  laboratory 
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Fig.  6. —  Laboratory  spectrum  of  hydrogenated  amorphous  carbon  healed  to 
800  K  (Scott  1 997).  The  two  portions  of  the  spectrum  separated  by  the  vertical 
hatched  line  at  7.5  pm  were  obtained  separately  and  have  been  arbitrarily  nor¬ 
malized  here  to  the  same  vertical  range.  The  apparent  relative  strengths  of  the  two 
halves  arc  not  meaningful.  The  positions  of  the  bands,  however,  are  a  good  in¬ 
dication  of  how  the  wavelengths  of  the  various  modes  responsible  for  the  observed 
PAH  features  may  shift  as  aliphatic  bonds  are  attached  to  the  aromatic  hydrocarbons 
that  form  the  basis  of  free  PAH  molecules.  Specifically,  the  two  strongest  features 
at  6.32  and  U.45  /xm  have  shifted  from  6.23  and  1 1 ,26  /mi  in  more  normal  PAH 
spectra.  The  vertical  lines  mark  the  positions  of  the  6.85  and  7,25  /on  features  seen 
in  spectra  of  cooler  HAC  samples,  [See  the  electronic  edition  of  the  Journal far 
a  color  version  of  this  figure .] 

HAC  spectrum  is  consistent  with  the  wavelength  shifts  we  have 
observed  in  class  C  PAH  sources* 

3.5*  PAHs.  HAC,  and  the  Class  C  Spectrum 

Bonds  in  aromatic  and  aliphatic  species  have  resonances  at 
different  wavelengths  (see  Bellamy  1978  and  the  many  refer¬ 
ences  therein).  In  the  case  of  the  C-H  stretching  fundamental, 
the  shift  from  aromatic  to  aliphatic  is  always  to  the  red,  from  3,29 
to  3,40  /xm,  but  for  other  vibrational  modes,  the  shifts  can  be  in 
either  direction,  or  both,  depending  on  the  molecules  in  question. 
Wc  hypothesize  that  at  longer  wavelengths,  the  substitution  of 
aliphatic  sidegroups  on  the  edges  of  large  PAH  molecules  shifts 
the  out-of-plane  aromatic  bending  modes  to  the  red.  While  lab¬ 
oratory  data  of  large  PAHs  with  aliphatic  substitutions  arc  lack¬ 
ing,  data  for  some  of  the  simplest  PAHs  are  available,  and  as 
Figure  7  shows,  are  consistent  with  our  hypothesis* 

Using  data  from  the  NIST  Chemistry  WebBook,  maintained 
by  the  National  Institute  for  Standards  and  Technology,3  0  we  show 
in  Figure  7  die  spectra  of  naphthalene  and  1 ,2-dimethylbcnzene, 
which  is  benzene  with  two  adjacent  substitutions  of  methyl  side- 
groups.  The  two  spectra  in  Figure  7  were  obtained  from  vapor 
samples  using  a  Fourier-transform  infrared  spectrometer  at  a  res¬ 
olution  of  4  cm"1  {Coblcntz  Society,  Inc,  2005)*  The  napthalene 
was  at  a  temperature  of  518  K,  while  the  1  f2 -dimethyl benzene 
was  at  423  K,  Since  both  molecules  have  four  adjacent  H  atoms 
on  each  ring,  the  only  out-of-plane  C-H  bending  mode  they 
show  in  the  1 1  - 14  jxm  region  is  the  quartet  mode.  This  mode  has 
shifted  from  ) 2.8  pm  in  naphthalene  to  13.5  pm  in  the  simpler 
benzene  molecule  with  the  methyl  sidegroups,  consistent  with 
the  anticipated  effect  of  the  substitution  of  a  methyl  sidegroup. 
However,  this  spectral  difference  may  also  result  from  the  ab¬ 
sence  of  the  second  ring  in  naphthalene,  making  further  analysis 
of  larger  molecules  highly  desirable. 

10 
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Fiti.  7.- —  Laboratory  spectra  of  two  simple  PAHs:  naphthalene  and  benzene 
with  two  adjacent  methyl  side  groups  (  from  Coblentz  Society,  Inc,  2005).  The 
substitution  of  aliphatic  si  degroups  on  benzene  adds  a  second  OH  stretching 
mode  at  3.4  /on,  produces  a  strong  absorption  at  6,8  /im,  shifts  the  quartet  out-of- 
plane  C-H  bending  mode  from  12,8  to  13.5  /im,  and  produces  more  spectral 
structure  in  the  vicinity  of  the  observed  6.85  and  7.25  /tm  aliphatic  features  (solid 
vertical  lines),  [fe  the  electronic  edition  of the  Journal for  a  color  i  version  of this 

Mure) 

We  also  note  that  the  presence  of  aliphatic  C-H  stretching 
modes  in  the  methyl  sidegroups  on  the  benzene  molecule  has 
added  emission  at  3.4  //nt  not  present  in  the  naphthalene.  Fur¬ 
thermore,  in  the  6-8  pm  region,  the  strongest  emission  feature  in 
naphthalene  is  at  6,6  /im,  but  a  stronger  feature  at  6.8  fim  is 
apparent  in  the  methyl -substituted  benzene.  As  described  by 
Bellamy  (1978),  the  carbon  skeletal  breathing  modes  in  aromatics 
cluster  in  the  6.2  to  6.9  /on  range,  while  similar  modes  from 
methyl  (CH3)  and  methylene  (CH2)  groups  cluster  in  the  6.8- 
15  fim  region.  As  noted  in  §  2.4,  we  may  have  detected  these 
methyl  and  methylene  features  in  emission  in  some  of  the  class 
C  PAH  spectra,  including  HD  1 00764. 

Laboratory  HAC  spectra  by  Furton  et  ah  { 1999)  dearly  show 
emission  features  at  6.85  and  7.25  //in.  The  HAC  spectrum  in 
Figure  6  only  shows  a  single  emission  feature  at  7.0  /jhi  at  about 
half  the  strength  of  the  6.3  /im  feature.  The  differences  in  the 
spectra  probably  arise  in  how  the  samples  were  treated  and  mea¬ 
sured.  Furton  et  al  ( 1999)  obtained  absorption  spectra  from  sam¬ 
ples  at  300  K,  while  the  HAC  spectrum  in  Figure  6  was  observed 
in  emission  at  —800  K.  This  higher  temperature  has  probably 
modified  the  simple  groups  responsible  for  the  6,85  and  7.25  /im 
features,  and  it  has  also  enhanced  the  C-C  modes  at  6.2  /im, 

HAC  samples  have  a  number  of  degrees  of  freedom,  includ¬ 
ing  the  ratio  of  aliphatic  and  aromatic  bonds  and  the  size  distri¬ 
bution  and  structure  of  the  embedded  PAHs,  among  others,  so  it 
is  not  surprising  that  no  HAC  sample  matches  the  observed  as¬ 
trophysical  spectra  in  detail.  It  is  worth  recalling  that  no  labo¬ 
ratory  PAH  spectrum  provides  a  precise  match  to  astrophysical 
PAH  spectra  either.  Available  data  arc  consistent  with  the  hy¬ 
pothesis  that  increasing  the  ratio  of  aliphatic  to  aromatic  bonds  in 
a  hydrocarbon  mixture  explains  the  differences  between  class  B 
and  class  C  PAH  spectra. 

We  propose  that  complex  astrophysical  hydrocarbons  are 
originally  synthesized  as  some  form  of  HAC,  When  this  material 
is  exposed  to  harsh  radiation  fields,  the  aliphatic  bonds  are  de¬ 
stroyed  and  only  free  PAHs  remain.  Scott  et  ah  ( 1997b)  pro¬ 
vide  evidence  that  this  process  could  occur.  They  exposed  HAC 
mixtures  to  UV  radiation  and  found  that  the  dissolution  of  the 
HAC  led  to  the  production  of  large  PAH  molecules  with  attached 


aliphatic  sidegroups.  Hu  et  ah  (2006)  have  recently  shown  that 
these  mixtures  of  aromatic  and  aliphatic  bonds  have  the  poten¬ 
tial  for  reproducing  the  observed  features  in  the  6-9  /mi  range. 
These  modified  PAHs  have  less  than  a  few  hundred  atoms  and 
still  respond  to  the  absorption  of  single  photons  like  more  typical 
free  PAHs  do. 

Two  separate  processes  are  responsible  for  the  PAH  spectra 
we  observe.  First,  the  PAHs  must  be  freed  from  the  aliphatic 
matrix  in  which  they  are  initially  embedded.  Second,  the  features 
must  be  excited.  This  latter  step  requires  less  energy  than  that 
needed  to  completely  strip  the  PAHs  of  all  aliphatic  side-chains. 
The  positions  of  the  observed  PAH  features  depend  on  how 
completely  the  PAHs  have  been  freed  from  the  underlying  ma¬ 
trix.  In  the  class  C  PAH  sources,  the  aliphatic  side  groups  have 
only  been  partially  destroyed. 

The  sample  of  isolated  HAeBc  stars  represents  a  stage  in  the 
evolution  of  PAHs  between  class  C  and  class  B.  The  7.85  /mi  fea¬ 
ture  is  shifted  measurably  to  7.9- 8.0  /im  or  beyond  (Sloan  et  al. 
2005;  Keller  et  al.  2007).  Here  again,  the  PAHs  may  not  be  en¬ 
tirely  stripped  from  the  underlying  aliphatic  matrix,  since  the  A 
and  F  stars  exciting  the  PAH  spectrum  may  not  produce  enough 
phoLons  with  sufficient  energy  to  break  all  of  the  aliphatic  bonds, 

3.6.  The  Problem  of  the  Reflection  Nebulae 

Reflection  nebulae  typically  do  not  show  class  C  PAH  spec¬ 
tra,  even  though  the  exciting  source  can  be  as  cool  as  6800  K 
(Uchida  et  al.  1998).  The  two  reflection  nebulae  in  the  sample 
examined  by  Peeters  et  al.  (2002)  are  both  class  A,  as  are  the 
three  spectra  presented  by  Uchida  et  al.  (2000).  The  objects 
studied  by  Bregman  &  Tcmi  (2005)  include  one  from  Uchida 
et  al.,  and  all  three  showed  variations  in  the  position  of  the 
7, 7- 7. 9  pm  feature  with  radiation  field,  but  not  enough  to  move 
them  even  as  far  as  a  true  class  B  spectrum. 

Wc  suggest  that  while  the  PAHs  in  these  complex,  evolving 
environments  are  now  subject  only  to  cool  radiation  fields,  they 
may  have  been  exposed  in  the  recent  past  to  a  large  dose  of  more 
intense  ultraviolet  radiation.  Reflection  nebulae  are  typically  dif¬ 
fuse  and  are  associated  with  complicated  star- forming  environ¬ 
ments.  This  diffuse  nature  means  that  the  PAHs  arc  poorly  shielded 
from  both  current  and  past  radiation  fields.  Furthermore,  PAHs  in 
a  diffuse  environment  arc  not  able  to  regenerate  any  broken  ali¬ 
phatic  bands. 

The  typical  PAH  emission  source  is  a  PDR  in  a  star  formation 
region,  and  the  emitting  materials  in  these  environments  must 
have  a  complicated  history  of  UV  irradiation.  This  history  blurs 
the  link  between  the  current  PAH  composition  and  the  current 
radiation  field,  which  is  probably  why  the  correlation  illustrated 
m  Figure  5  has  not  been  noticed  before.  Only  in  more  isolated 
sources,  such  as  those  examined  here,  is  the  link  clear. 

4.  CONCLUSIONS 

The  PAH  spectrum  of  HD  100764  adds  to  the  small,  but  grow¬ 
ing  sample  of  class  C  PAH  spectra,  as  defined  by  Peeters  et  al, 
(2002),  It  has  one  of  the  coolest  exciting  stars  and  some  of  the 
most  redshifted  features  in  the  class.  Since  all  seven  of  the  known 
class  C  spectra  arc  excited  by  relatively  cool  stars  of  spectral 
class  F  or  later,  we  argue  that  the  hydrocarbons  in  these  systems 
arc  unusual  because  they  have  not  been  exposed  to  much  ultra¬ 
violet  radiation. 

We  hypothesize  that  carbonaceous  materials  are  synthesized 
as  large  HAC  conglomerates  and  that  aliphatic  bonds  are  sub¬ 
sequently  broken  in  harsh  radiation  fields.  We  suggest  that  PAH 
spectra  evolve  over  time;  the  hydrocarbons  in  class  C  objects  are 
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relatively  protected  and  unprocessed,  while  class  B  and  A  PAHs 
have  been  exposed  to  more  energetic  photons  and  are  more 
processed. 
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